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The resonant activation phenomenon (RAP) in a discrete system is studied using the master 
equation formalism. We show that the RAP corresponds to a non-monotonic behavior of the fre- 
quency dependent first passage time probability density function [pdf). An analytical expression for 
the resonant frequency is introduced, which, together with numerical results, helps understand the 
RAP behavior in the space spanned by the transition rates for the case of refiecting and absorbing 
boundary conditions. The limited range of system parameters for which the RAP occurs is discussed. 
We show that a minimum and a maximum in the mean first passage time (MFPT) can be ob- 
tained when both boundaries are absorbing. Relationships to some biological systems are suggested. 

PACS number: 05.40. Ca, 82.20.Db, 02.50.Cw 



1. INTRODUCTION 

Noise induced escape of a particle from a potential well 
has been a fundamental way for describing various pro- 
cesses in biology, chemistry, and physics, since the sem- 
inal work of Kramers ll. More recently, it has been 
suggested that for some systems the potential itself fluc- 
tuates in time. A few examples are: the transport of ions 
and bio-polymers through membrane channels en- 
zymatic kinetics , and the rebinding of ligand-protein 
complexes 0, 0| • The basic formulation of noise induced 
escape from a fluctuating environment is obtained by 
making the potential term of the stochastic differential 
equation (e.g. a white noise overdamped Langevin equa- 
tion) change with a frequency 7 between two states. Do- 
ering and Gadoua showed that for a fluctuating system, 
the MFPT, T, from a reflecting boundary to an absorbin 
boundary, may show a minimum as a function of 7 
The occurrence of a minimum in t(7) was termed the 
resonant activation phenomenon (RAP). This has been 
followed by an extensive theoretical work to understand 
the nature of the RAP 9]- 16], along with experimental 
efforts to find systems displaying RAP 0, 0| . The the- 
oretical works have been mainly focused on checking the 
effect of different potentials on the RAP. 

In this paper we study the discrete case RAP using 
coupled master equations (ME). We show that the RAP 
is only one of the properties that stem from the non- 
monotonic behavior of the frequency dependent first pas- 
sage times (FPT) pdj, Fj{t), and which are related to 
frequency dependent minima in the first and higher mo- 
ments of F^{t). We introduce an analytical expression 
for the dependence of the frequency that minimizes t(7) 
on the system transitions rates. We show that the RAP 
is obtained only when certain conditions imposed on the 
transition rates are fulfilled. Analyzing these conditions 
we come up with an instructive understanding regard- 
ing the nature of the RAP. In addition, a novel behavior 
of the MFPT is obtained when changing the reflecting 





boundary into an absorbing one: the coexistence of a 
minimum and a maximum in t(7). 

2. THEORETICAL FRAMEWORK 

We describe the escape from a system that fluctuates 
between two configurations A and B, using the coupled 
ME: 

A — 17 I7 
I7 B I7 

(1) 

— ? — ? 

P A{t) {P sit)) is an n-dimensional column vector, whose 
j element is the pdf to occupy site j of the A (B) con- 
figuration at time t. The transition between each site 
j in one configuration and its counterpart in the second 
configuration occurs with a flipping frequency 7, see Fig. 
1. I is the unit matrix of n dimensions introduced in Eq. 
(1) to indicate the configurational coupling. Movement 
along each of the configurations, A and B, is governed by 
the square n-dimensional tridiagonal propagation matri- 
ces A and B, respectively, whose elements are the transi- 
tion rates (Fig. 1). The choice of the matrices A and B 
corresponds to an equivalent choice of potential profiles 
and boundary conditions in the continuum case. In what 
follows we set a reflecting boundary at site j—n, and an 
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Figure 1: Schematic illustration of the exit problem in a fluc- 
tuating environment for n — 2. For an invariant system 
ai,o=a2,i=a/, ai,2=t[6, and &i,o=&2,i=fe/, 6i,2=&b. For the 
particular invariant birth-death system ab=0, and 6/=0. 
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absorbing boundary, as a trap, at site j=0, unless other- 
wise indicated. Fig. 1 shows a schematic illustration of 
the coupled system for 71 = 2. 

The FPT pdf is defined by F^{t) = 9(1 - S^(t))/dt, 
where S^it) is the survival probability; namely, the prob- 
ability of not reaching the site j=0 until time t. Sj{t) 
is obtained by summing the elements of the vector that 
solves Eq. (1), S^{t) = [/2„Ee°*E-i P2n(0). Here U2n 
is the summation row vector of 2n dimensions, P2n(fi) 
is the initial condition column vector, [-P2n(0)]j = 
(Sr^ j Pji Q+Sx+n,j Pb ,0) : where x is the initial site, and the 
process starts in the A (B) configuration with probability 
Pa,o {Pb,o)- Unless otherwise specified, we use x=n as 
a starting site, and Pa,o=-Pb,o=1/2, as suggested from 
the single configurational flipping frequency. The defi- 
nite negative real part eigenvalues matrix D is obtained 
through the similarity transformation: D = E^^HE, 
where H is the matrix given on the right hand side of 
Eq. (1), and E and E^^ are the eigenvectors matrix, and 
its inverse, of H. 

3. RESULTS AND DISCUSSION 

We start by computing Fj{t) for an invariant birth- 
death system. By an invariant system we mean that 
the transition rates are independent of the site index, j; 
namely, mj,j_/+i=m / /b for configuration M, where 
represents and bf/t transition rates, and M stands 
for A and B. A birth-death system means that the move- 
ment in each configuration occurs only in one direction, 
i.e. ab=bf=0. Clearly, the term birth-death indicates 
that the particle [when simulating Eq. (1)] can move 
only towards its "death" (the trap) when it is subjected 
to the dynamics of the A configuration, and in this sense, 
when flipping to the B configuration occurs it is "born" 
(or "resurrected"). Therefore, for the birth-death sys- 
tem, the fluctuations are between a configuration which 
acts as a "barrier", the birth configuration, and a con- 
figuration acting as a "valley", the death configuration. 
Note that a single-rate (namely, af=bb) invariant birth- 
death system is similar to the system studied by Doering 
and Gadoua , where the derivatives of the two linear 
potentials are sign opposite and equal in the absolute 
values. 

Fig. 2 shows Fj{t) for a single-rate invariant birth- 
death system and n=10. At short to intermediate times, 
Fj{t) displays a peak that shifts towards larger times 
with 7. This peak represents the exiting of the initial 
population of configuration A, 6x,jPA,Q, for small 7, and 
the overall initial condition for large 7. At longer times 
and intermediate 7, a minimum in F^(t) appears as a 
function of 7 that represents the fastest exit mainly of 
the initial B population. Accordingly, the minimum in 
the MFPT t{j) is a consequence of the shape of Fy{t), 
and is therefore reflected in higher moments of Fry{t) as 
well. 
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Figure 2: F.y(t) as a function of t and 7 on a log- log-log plot 
(natural log), for a single rate invariant birth-death system 
with n=10, and a/=6i,=l. Also shown is F^{t) projections. 

To study the RAP we start by computing the 
MFPT. In general, the s*'' moment of F{t) is ob- 
tained by inverting matrix H: = t''F{t)dt = 

s!C^2n(— H)^^P2n(0). U siug the projection operator 
techniques for H^^ blocks r reads 

t = 7C/„(Ca + Cb)P„(0)- 

-C/„(CaA + CbB)P„(0)/2, (2) 

where Ca=[AB-7(A + B)]-i, Cb = [BA-7(A + B)]-i, 

and [P n{Q)]j=Sx,j- From Eq. (2) one can calculate 
t(7) for the limiting cases 7^0 and 7 — > 00. For 
7 — > 0, r is the average of the MFPT of the uncou- 
pled configurations, A and B, t={ta + tb)/'2, where 
TM = -c7m-1P„(0) is the MFPT of configuration M. 
For 7 ^ 00, r is the MFPT of an averaged fully cou- 
pled system, namely, r=— U „ ( ) ^ ^ P „ (0) . These are 
the expected hmiting behaviors of the MFPT 
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RAP is expected for intermediate hipping fre- 
quencies. 

To obtain an analytical expression for the frequency 
that minimizes t(7), 7mm, we search for an extremum 
point (a minimum) of the function r(7), for an invariant 
system and n=2. We find jmin to be a sum of two terms: 



7mm = hmin,! > 0] + [7mm, 2 > 0], 



(3) 



where the notations on the right hand side of Eq. 3 mean 
that each of the terms must be non negative to contribute 



to 7„ 



and 



„, and 

af{af ^2bb) ~bf{bf -2ab) 
- a/^b) - bf{a} - bfttb) 



af{af + 26fo) — bf{bf + 2a^,) 



(5) 
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We note that the smallest system that exhibits the RAP 
requires a three site system, which is a specific case of 
the system shown in Fig. 1, with, for example, 61^2 = 
62.1 00. However, in what follows we consider systems 
with finite transition rates. 

For the birth-death system Eq. (3) reduces to: 



and 7=(rA) ^ into t(7), results in r(7 = "/.mm) ~ 2/n2, 



7" 



2 - Uf/bb' 



(6) 



The simple form of Eq. (6) provides an insight into 
the nature of the RAP. It immediately implies the re- 
quirement af/bf, < 2 for RAP to occur. For bb ^ a/, 
Jmin={TA)^^ , where ta is the first moment, s—1, of 
F{t) for a death system, rj^ = {n)s/a'p where {n)s = 
(n + s — l)!/(n — 1)!. This optimal frequency means that 
Sx.jPas) has exited the interval, on average, while the 
first configurational transition occurred, and the same 
holds for SxjPsfi, for the second configurational tran- 
sition. Because the probability (particles) can exit the 
interval only when it is subject to the A configuration dy- 
namics, a situation where the A configuration is empty 
but not the B configuration, means a "waste" of time 
with regards to fastest interval exiting. This is the case 
for 7 < jrnin- For 7 > jrnin, not all 5x,jPA,o exited the in- 
terval, while the first configurational transition occurred, 
meaning that another cycle of flipping is required to exit 
the system. This leads again to a "waste" of time with 
regards of fastest interval exiting. At 7 = {ta)~^, only 
one configurational change occurs, and costs the minimal 
time for exiting the interval. 

The special feature that Uf/bb < 2 is needed for RAP 
suggests that the rate along the birth configuration must 
be, at least, as fast as those along the death configura- 
tion for the RAP to be obtained. For an invariant birth- 
death system to show the RAP, the ratio af/bb must fulfil 
df/bb < 3 asymptotically, which is demonstrated in Fig. 
3. Note that Fig. 3 spans both degrees of freedom of 
the invariant birth-death case, the size n, and the ratio 
af/bb- Scahng the time, t = tbb, leads to dimensionless 
rates af/bb and j/bb- 

We emphasize that even for the simple invariant birth- 
death system that fulfills the demand that the B config- 
uration rate is much larger than of configuration A rates, 
the relation -ymin = (t'a)"^ might not be satisfied. To see 
that we calculate -ymin for a general birth-death system 
and n=2 (see Fig. 2): 



7n 



2 - ^01,002,1/61,2 



(7) 



01.002.1/(0.1.0 

-1 



which for ^ai,oa2,i/J'i,2 
002,1/2, where (ta)~ 
over, when ai,o » 02,1, (ta) 
changed, and is much larger than (ta) 
more than one configurational change occurs at the op- 
timal flipping frequency. Reverse substitution "f=jmin 



0, reduces to 7mm = 
02,1 )• More- 
02,1, 7mm is un- 
^, implying that 



and 



0.8. Note that for both values of 



t(7 = (ta) ^) 

7, t{j) is independent of 61,2 to be compared with 
r(7 ^ 00) w s:^^^- Imin has, therefore, a general 
functional form that not necessarily coincides with the 
MFPT of the faster configuration. 

Going beyond the birth-death system, we first consider 
a case for which B=AA. From Eqs.(3)-(5) we have, for 
n=2, "imin < for any positive A. From numerical cal- 
culation we find that for n > 2 there is no real positive 
"imin- Both the analytical and the numerical results im- 
ply that a system for which A and B commute does not 
exhibit the RAP. 

The next case to be checked for the occurrence of the 
RAP is obtained by setting, af=bb=k. Using t = kt, the 
system transition rates are dimensionless, measured in 
units of k. This procedure leads to the reduced Eqs.(4- 
5): 



'ymin,l/ k '■ 



f2[3- 



3u 



)[2u-v\ - 1 



and 



^minfi/k — 



2,-v[2u + v\ 



(8) 



(9) 



where u=ab/k and v—bf/k. Fig. 4a shows "fmin/k as a 
function of v for u=l/2. "fmin,i/k displays a maximum 
at Vmax, which is easily recovered from Eq. (8). For 
v < Vmax, IminA/k iucrcascs, which reflects the increase 
in the relative ability of the B configuration to "help" 
the fastest exiting of configuration A. On the other hand, 
the decrease in "fmin,i/k for v > Vmax, implies that the B 
configuration movement towards the absorbing end be- 




Figure 3: "fminTA/bb as a function of n and the ratio af/bb- 
When ttf/bb > 3 the RAP does not exist, for any n. 
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Figure 4: a: ymin/k for af=bi,—k, and u=l/2, as a function 
of «. b: 'ymin.i/k as a function of < ?; < 1. c: 'ymin/k for 
11=1 as a function of v. d: 'ymin/k for it=5 as a function of ii. 



comes fast enough "to stand on its own" for the accom- 
pUshment of this task. A resonant-free zone occurs in the 
region where ta ~ tb, and is followed by a short reso- 
nant region, where both configurations are trap oriented, 
namely, u < 1 and v > 1. Fig. 4b shows that in the range 
< u < 1 ^min.i/k is non-monotonic. Fig. 4c shows for 
u—1 the only non-zero 7mm,i/fc > which diverges as 
when u — !■ 1, because for these system parameters 
B=A. For u > 1 and v < 1, ^min/k has two asymptotic 
lines, which define a resonant-free region. This happens 
when B « AA. For u ^ v > 1, jmm/k ~ v/2. These 
two features are demonstrated in Fig. 4d. 

The above analysis is of importance since the cou- 
pled ME with n=2 can be used to model the kinetics 
of a conformationally changing enzyme. Such extended 
Michaelis-Menten models are appropriate for describing 
experiments performed on a single molecule level [5]. If 
we assume that there are two enzyme conformations, a 
specific stage of the enzymatic activity can be described 
by Fig. 1. From Eqs. (8) and (9), and more generally 
Eqs. (4) and (5), a relationship between the reaction 
rates and the conformational flipping frequency can be 
established for an optimal enzymatic activity. In addi- 
tion, changes in the flipping rate value near the resonant 
frequency, which can be achieved for example by binding 
of other molecules to the enzyme, provide a simple and 
efficient mechanism for regulating the enzymatic activity, 
which is a well known issue in biology 

Finally, we study a system for which both ends are 
absorbing; namely, the reflecting boundary is replaced 
by an absorbing one, and the escape process starts at 
the middle site {n—7, and x=4). The coupled invariant 
conflgurations are taken to have an opposite bias, i. e. 
the transition rates of the B configuration are 6/=0.175, 
bb=l, which give rise to a left side bias as defined in Fig. 



1, whereas in the A configuration a/=3.5, and Ub < af 
give rise to a right side bias. A global minimum and a 
global maximum in t{'-^)/t{-^ — > oo) can occur (Fig. 5b). 
The minimum and maximum appear at the neighborhood 
of the points jminTB = 1, and jmaxTA = 1, respectively. 
This behavior is sensitive to the value of a?, (Fig. 5a). 
For a;, — > the global extremum points reduce to local 
extremum points. When a{, — > o/, t(7)/t(7 —^ oo) is a 
monotonically increasing function of jta- 

We note that these boundary conditions for a fluctuat- 
ing system have been used to describe the translocation 
of a single stranded DNA through a conformationally 
changing nanopore fj]. For this case no resonance oc- 
curred because of the physical conditions that imposed 
the relation B=AA. However, for systems that are de- 
scribed by matrices A and B that do not commute and 
for two absorbing ends, a change in 7 in the vicinity 
of the extremal points, leads to a drastic change in the 
average time during which the system is occupied, and, 
therefore, emphasizes the importance of the frequency of 
fluctuation as a control parameter. 



4. CONCLUSIONS 



To conclude, in this paper we revisited the resonant ac- 
tivation phenomenon. We studied the origin of the RAP 
and the requirements under which this phenomenon can 
be observed. We showed that for a single rate invariant 
birth-death system the RAP is a consequence of a general 
phenomenon, which is a non-monotonic behavior of (t) 
along the frequency axis for large time. We characterized 
the conditions for which an invariant birth-death system 
exhibits the RAP, and broaden these conditions by ex- 
amining more general systems. Relationship between the 
RAP and biological activity was suggested. In addition, 
we introduced a new property of the MFPT, the coex- 
istence of a minimum and a maximum in the flipping 
frequency dependent MFPT, r(7). 




Figure 5; a: A log-log-log plot (natural log) of t"(7)/t(7 — > 
00) as a function of ai, and 7rA , for two absorbing ends in 
an invariant system, with a/=3.5, 6/=0.175, 6i,=l, n=7, and 
x=4. b: Profile of the left figure for ai,=0.25 exposes a global 
minimum and a global maximum in t(7)/t(7 — > 00). 
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